import pandas as pd
import numpy as np

def add_noise_and_save(file_path, output_path, sigma=0.0042):
    """
    给 t90 添加正态扰动，并保存新数据
    """
    # 读取原始数据
    df = pd.read_excel(file_path)

    # 添加扰动列
    df['Y_noisy'] = df['Y染色体浓度'] + np.random.normal(loc=0, scale=sigma, size=len(df))

    # 保存为新的文件
    df.to_excel(output_path, index=False)
    print(f"已保存扰动后的数据: {output_path}")

if __name__ == '__main__':
    add_noise_and_save('data.xlsx', 'data_noisy.xlsx', sigma=0.0042)
